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Abstract 



A numerical study of low- lying glueball masses of compact U(l) lattice gauge 
theory in (2+1) dimensions is performed using Standard Path integral Monte 
Carlo techniques. The masses are extracted, at fixed (low) temperature, from 
simulations on anisotropic lattices, with temporal lattice spacing much smaller 
than the spatial ones. Convincing evidence of the scaling behaviour in the an- 
tisymmetric mass gap is observed over the range 1.4 < /3 < 2.25. The observed 
behaviour is very consistent with asymptotic form predicted by Gopfert and 
Mack. Extrapolations are made to the "Hamiltonian" limit, and the results 
are compared with previous estimates obtained by many other Hamiltonian 
studies. 
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I. INTRODUCTION 



Compact U(l) lattice gauge theory in three- dimensions is the most well-understood non- 
trivial gauge theory. The theory exhibits some important similarities with QCD [16], such 
as confinement [2,3] and chiral symmetry breaking [4]. Because of its simplicity, it has 
been a good testing ground for new methods and algorithmic approaches. While the Eu- 
clidean Path Integral Monte Carlo simulations [5] have been very successful in the study of 
non-perturbative lattice gauge theories and is currently the preffered technique for studying 
QCD at low energies, there there has been some progress in the development of analytic and 
numerical approaches in the Hamiltonian formalism. The Hamiltonian version of the model 
has been studied by many methods: some recent studies include series expansions [6] , finite- 
lattice techniques [7], the t-expansion [8,9], and coupled-cluster techniques [10-12], as well 
as Quantum Monte Carlo methods [13-19]. Quite accurate estimates have been obtained 
for the string tension and the mass gaps, which can be used as comparison for our present 
results. The finite-size scaling properties of the model can be predicted using an effective 
Lagrangian approach combined with a weak-coupling expansion [20], and the predictions 
agree very well with finite-lattice data [16]. Here our aim is to use standard Euclidean path 
integral Monte Carlo techniques for an anisotropic lattice, and see whether useful results 
can be obtained in the anisotropic or Hamiltonian limit. 

II. COMPACT U(l) MODEL IN 3 DIMENSIONS 

The anisotropic Wilson gauge action for compact U(l) model in (2-1-1) dimensions has 
the form [21] 

S = (isY.P^3{r)+fitY.Prtir) (1) 

r,i<j r,i 

Where and Pa are the spatial and temporal plaquette variables respectively. In the 
classical limit 

^ 0, % 9 
1 11 

A = ^ = -7^ (2) 
e^at At 

where (3 — l/g^ {g"^ — ae^ in (2-|-l)D)) and At = at/ag is the anisotropy parameter, is the 
lattice spacing in the space direction, and at is the temporal spacing. In the weak-coupling 
approximation, the above action can be written as 



S = /3 



E E ( 1 - ) + ^ E ( 1 - ^^t(r) 



r i<j 



(3) 



In the limit At — > 0, the time variable becomes continuous and one obtains the Hamiltonian 
hmit of the model (modulo a Wick rotation back to Minkowski space) . 

Among other features, antisymmetric mass gap is a fundamental quantity of interest in 
U(l) model. On an isotropic lattice, the model reduces to the simple continuum theory of 
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non-interacting photons in the naive continuum hmit at a fixed energy scale [22] but if one 
renormahzes or rescales in the standard way so as to maintain the mass gap constant, then 
one obtains a confining theory of free massive bosons. Gopfert and Mack [23] proved that 
in the continuum hmit the theory converges to a scalar free field theory of massive bosons. 
They found that in that limit the mass gap behaves as 




arno = W ^ exp( -'f;(0)) (4) 



where v{0) = 0.2527 is the Coulomb potential at zero distance for the isotropic case. 

The behaviour of the mass gap in the anisotropic case will be similar to equation (4) 
Generalizing discussions by Banks et al [24] and Ben-Menahem [25], we find that the expo- 
nential factor takes exactly the same form in the anisotropic case. The only difference is 
that the lattice Coulomb potential at zero spacing for general At is 

_ r At , , 

'"^^^ " (27r)3 4[sin2(A:o/2) + AT^{sm\ki/2) + sm\k2/2))] ^' 

_ j 0.2527 (At = 1) . . 

\ 0.3214 (At = 0) 

The above result is based on dilute gas approximation, which is justified in the Euclidean 
case, but not in the Hamiltonian limit [24,25]. 

The expected finite-size scaling behaviour of the mass gap near the continuum critical 
point in this model is not known; but Weigel and Janke [26] have performed a Monte 
Carlo simulation for an 0(2) spin model in three dimensions which should lie in the same 
universality class, obtaining 

M ~ 1.3218/L (7) 

for the magnetic gap. 



III. MONTE CALRO SIMULATIONS AND GLUEBALL MASSES 

Using the anisotropic Wilson action, we perform Monte Carlo simulations on a finite 
lattice of size N'^ x iV^, where Ng is the number of lattice sites in the space direction and 
Nr in the temporal direction, with spacing ratio Ar = at/ as- By varying Ar it is possible 
to change at, while keeping the spacing in the spatial direction fixed. The details of the 
algorithm for updating are given elsewhere [27]. The simulations were performed on lattices 
with Ng — 16 sites in each of the two spatial directions and Nt — 16 — 64 in the temporal 
direction for a range of couplings P — 1 — 3. 



IV. RESULTS AND DISCUSSION 

Figure 1 shows the estimates for the glueball masses of 0++ and channels against (3 for 
the Euclidean case At = 1.0. As expected the masses are seen to vanish exponentially with 
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/3 as we approach the continuum hmit (high f3). To show the evidence of scahng behaviour 
of antisymmetric glueball masses , we used the predicted asymptotic form, equation 4, with 
an additional normahzation constant; 



where ci = 5.23 ±0.11 when adjusted to fit the data. Figure 2 shows the predicted scahng 
behaviour. The sohd hne on the graph is a fit to the data over the range lA < (3 < 2.25. 
The slope, 2.48 ± 0.09, of the data matches the predicted asymptotic form very nicely, but 
the coefficient is too large by a factor of 5.2. It would again be interesting to explore whether 
this discrepancy could be due to quantum corrections. The finite-size scaling behaviour is 
shown by the dashed line in Figure 2. It can be seen that the Euclidean mass gap should 
not be effected by finite-size corrections until /3 > 2.8. 

To check the consistency of our method, we plot the dimensionless ratio of the the 
antisymmetric mass gap over the square root of the string tension against (3 together with 
the results of Teper [29] in Figure 3. The agreement is excellent. The solid line gives 
the ratio of the fits in Figures 2 and Fig. (8) in [28], and shows how this ratio vanishes 
exponentially in the weak-coupling limit, whereas in four-dimensional confining theories it 
goes to a constant. 

Figure 4 shows the behaviour of the glueball masses as function of Ar^ for /3 = \/2. The 
extrapolation to the Hamiltonian limit is performed by a simple cubic fit in powers of Ar^. 
In the limit At 0, we reproduce the earlier series estimates of Hamer et al [6] for the 
0++ and states. The estimates of the antisymmetric mass gap, for various extracted 
from the extrapolation to the At — > hmit are graphed as a function of P in Figure 5. The 
solid line is a fit to the data for 1.4 < /5 < 2.25 of the form (4), with v{0) = 0.3214 and 
Ci = 5.50 ± 0.24, which is similar to the coefficient found in the Euclidean case. The fitting 
parameters, the slope, 3.10 ± 0.16 and the intercept, 3.61 ± 0.25, of the scaling curve are in 
excellent agreement with the results obtained from the other Hamiltonian studies and are 
hsted in table I. The dashed line represents the finite-size scaling behaviour, equation (7), 
which we assume holds in the Hamiltonian limit also, for want of better information. It 
can be seen that the finite-size corrections are predicted to dominate for f3 > 2.2, but the 
date are not accurate enough at weak couplings to establish whether this is really the case. 
Figure 6 shows the our Monte Carlo estimates for antisymmetric mass gap as a function of 
p. Also are shown the results from previous strong-coupling series extrapolations [6] and 
quantum Monte Carlo calculations [16]. It can be seen that our present results agree with 
previous estimates but are less accurate. 

Finally, Figure 7 displays the behaviour of the dimensionless mass ratio. 



for the Euclidean case At = 1. As in the (3-|-l)D confining theories, we may expect that 
quantities of this sort will approach their weak-coupling or continuum limits with corrections 

of 0(l/ae//), where a^ff is the effective lattice spacing in 'physical' units when the mass gap 
has been renormalized to a constant. Hence for our present purposes we define agff from 
equation (4) as 




(8) 



Rm = M(0++)/M(0~) 



(9) 




(10) 
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with v{0) = 0.2527 for the Euchdean case. The mass ratio is plotted against Oe// in Figure 
6. At weak couphng, we expect the theory to approach a theory of free bosons [23] so that 
the symmetric state will be composed of two bosons and the mass ratio should approach 
two. Our results show that as a^/f goes to zero, the mass ratio rises to around the expected 
value of 2.0. A linear fit to the data from 0.08 < Ce// < 0.35 gives an intercept However, we 
note that the last one of our estimates, together with two from Teper [29], lie considerably 
above Rm = 2. In the bluk systems, of course, the ratio cannot rise above 2, because it is 
always possible to construct a 0"*"+ state out of two mesons, been included in the fit. 
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TABLES 



Source 


Co 


Cl 


Villain (Hamiltonian) 


3.17 


2.18 


Lana [30] 


2.05 


2.49 


Hamer and Irving [31] 


2.65 


3.07 


Hamer, Oitmaa and Weihong [6] 


2.71 


3.13 


Dabrighaus, Ristig and Clark [10] 


2.40 


3.13 


Fang, Liu and Guo [11] 


2.50 


2.94 


Morningstar [9] 


2.61 


2.97 


Heys and Stump [32] 


2.48 


3.13 


Mcintosh and Hollenberg [17] 


2.50 


2.91 


Xiyan, Jinming and Shuohong [33] 


2.50 


2.95 


Darooneh and Modarres [34] 


2.20 ±0.03 


2.96 ± 0.05 


This work 


3.10 ±0.16 


3.61 ±0.25 



TABLE I. Comparison of the coefficients cq ( the negative of the slope of scaling curve) and ci 
(intercept of the scaling curve on the ln(M/ yfP) axis) in the weak coupling limit for antisymmetric 
mass gap. 



7 
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FIG. 1. The glueball mass estimates of as a function of /3. The estimates for 0++ and 
channels are shown by sohd triangles and open circles respectively. 
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FIG. 2. The scaling behaviour of the antisymmetric mass gap against /? at At = 1.0. The sohd 
hne is a fit of the form eq. (4) . The errors are smaller than the symbols. The dashed line shows 
the finite size scaling behaviour [26] 
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FIG. 3. The dimensionless ratio Mj\J~K as a function /?. Our estimates are shown by circles 
and solid triangles show the earlier results of Teper [29] . The solid curve represents the predicted 
weak-coupling behaviour. 
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AT" 

FIG. 4. Estimates of the masses of 0^"'' and glueballs against Ar^. Results at /? = \/2 for 
the 0^^ and are labeled by circles and triangle respectively. The solid and dashed curves are 
the cubic fits to the data extrapolated to the Hamiltonian limit. The series estimates of Hamer 
et al [6] in the limit Ar — 0, for symmetric and antisymmetric channels are shown as a star and 
diamond respectively. 
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FIG. 5. Hamiltonian estimates of the antisymmetric mass gap plotted as a function of (3. The 
solid curve is the fit to the data for 1.4 < /? < 2.25. The dashed line represents the finite size 
effects [26]. 
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FIG. 6. Hamiltonian estimates of the antisymmetric mass gap plotted as a function of /3. The 

solid curve represents the previous results from series expansion [6] and the dashed line represents 
the finite size effects [26]. The previous quantum Monte Carlo calculations [16] are shown as solid 
triangles. 
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FIG. 7. A graph showing estimates of the mass ratio Rmss a function of the effective spacing, 
tteff, at At = 1.0. Our present estimates are shown by the circles. The dashed hue is a hnear fit 
to the data over the range 0.08 < ag// < 0.35. The sohd triangles show the previous estimates of 
Teper [29]. 
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